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We consider a new type of cooling mechanism for a suspended nanowire acting as a weak link 
between two superconductive electrodes. By applying a bias voltage over the system, we show 
that the system can be viewed as a refrigerator for the nanomechanical vibrations, where energy 
is continuously transferred from the vibrational degrees of freedom to the extended quasiparticle 
states in the leads through the periodic modulation of the inter- Andreev level separation. The 
necessary coupling between the electronic and mechanical degrees of freedom responsible for this 
energy-transfer can be achieved both with an external magnetic or electrical field, and is shown to 
lead to an effective cooling of the vibrating nanowire. Using realistic parameters for a suspended 
nanowire in the form of a metallic carbon nanotube we analyze the evolution of the density matrix 
and demonstrate the possibility to cool the system down to a stationary vibron population of ~ 0.1. 
Furthermore, it is shown that the stationary occupancy of the vibrational modes of the nanowire 
can be directly probed from the DC current responsible for carrying away the absorbed energy from 
the vibrating nanowire. 

PACS numbers: 73.23.-b Electronic transport in mesoscopic systems; 

74.45.-l-c Proximity effects; Andreev reflection; SN and SNS junctions; 

85.85.-|-j Micro- and nano-electromechanical systems (MEMS/NEMS) and devices. 
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1. INTRODUCTION 

Nanoelectromechanical systems (NEMS) have over the 
last two decades been a very active field of both funda- 
mental and applied research. These systems have typi- 
cal dimensions on the nanoscale and combine electronic 
and mechanical degrees of freedom for novel applications. 
These include ultra sensitive mass detection [iHSl and 
position sensing [34^ using small mechanical resonators 
with a resolution that has been shown to approach the 
limit set by the Heisenberg uncertainty principle. Due to 
the small dimensions implied by NEMS they border on 
the world of quantum mechanics while still being macro- 
scopic in the sense that they can be fabricated litho- 
graphically. Nanoelectromechanical systems thus allow 
for the controlled fabrication of devices whose physics 
is ultimately governed by the fundamental limits set by 
quantum mechanics. Achieving control and understand- 
ing of these systems is therefore of great interest as doing 
so opens up a whole new toolbox for the design of high- 
performance apphcations 0I1|- 

Typically, nanoelectromechanical systems comprise a 
mechanical resonator in the form of a cantilever or a 
doubly-clamped beam coupled to an electronic system 
used for both actuation and detection. Due to the high 
mechanical vibration frequencies and exceptionally high 



quality factors recently achieved, these systems allow for 
very low energy dissipation and extreme sensitivity to ex- 
ternal stimuli. To make full use of their potential much 
research has recently focused on the possibility to effec- 
tively cool the mechanical subsystem in NEMS to its 
vibrational ground state, which would enable unprece- 
dented mass- and position-sensitivity of the resonators. 
Recently, O'Connell and colleagues [9jl showed that com- 
plete ground state cooling of a mechanical resonator can 
be achieved at very low temperatures. In Ref. Q the re- 
ported vibron occupancy factor is < 0.07, which means 
the resonator is in its ground state with a probability 
of 93%, well within the quantum regime. Here, we note 
however, that this effect is attributed to the high mechan- 
ical frequency of the resonator (in the GHz range) and 
the very low (25 mK) temperature of the dilution fridge. 
No active cooling of the resonator is needed as the low 
occupancy of its excited vibrational modes simply fol- 
lows from the Bose-Einstein distribution function at the 
equilibrium temperature. To date, the best experimental 
result for active cooling of a mechanical resonator lead 
to a vibron occupation factor of 3.8, as recently reported 
by Rocheleau et al. 10]. In Ref.tlCi] an external electro- 
magnetic field acts as an energy transducer for side-band 
cooling of the mechanical resonator (see also Ref. fl]|, 
where dynamical back-action is used to achieve similar 
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results) . 

In this paper we will consider a new experimental 
setup for ground state cooling of a mechanical resonator, 
a nanomechanical superconducting weak link as shown 
schematically in Fig. [TJ Below we show that this sys- 
tem can be viewed as a refrigerator for the nanomechan- 
ical vibrations, where the voltage-driven Andreev states 
serve as the refrigerant responsible for pumping energy 
from the nanomechanical vibrations to the thermostat 
of quasiparticle states in the leads. The possibility of 
such an energy transfer is based on the periodic, in time, 
modulation of the inter-Andreev level spacing and the 
coupling between these current carrying Andreev states 
and the mechanical degrees of freedom through an ap- 
plied magnetic field (we note that this coupling could also 
be achieved through an external electric field as outlined 
in Appendix E)) . At the start of each cooling cycle the 
Andreev level separation initially shrinks, which brings 
them into thermal contact with the vibronic degrees of 
freedom, enabling energy exchange between the two sys- 
tems. During this time, the work done by the bias voltage 
results in an "over-cooling" of the Andreev level popula- 
tion, which makes the heat transfer from the nanome- 
chanical to the electronic degrees of freedom possible. 
After a while the inter-Andreev level separation increases 
again, eventually reaching its initial value, at which point 
the Andreev states release their excess energy into the 
extended quasiparticle states of the leads. The cooling 
process suggested here thus models a nanomechanical re- 
frigerator where heat is continuously transferred along 
the chain; nanomechanical vibrations Andreev levels 
— >■ quasiparticle states, resulting in an overall effective 
cooling of the nanomechanical vibrations. 

Unlike many recent theoretical models in which elec- 
tromechanically induced ground state cooling of mechan- 
ical resonators have been reported [l^ - [l6j . the cooling 
mechanism suggested here relies neither on the above- 
mentioned side-band effect nor on dynamical back-action. 
Rather, the cooling is an inherent property of the de- 
vice which makes the suggested mechanism a promising 
candidate for truly quantum mechanical manipulation of 
mechanical resonators. 

The electromechanical coupling required for our cool- 
ing mechanism to work will be assumed to be due 
to an external magnetic field applied perpendicular to 
the length of the nanowire [13] ■ Displacements of the 
nanowire will then couple the mechanical and electronic 
degrees of freedom through the Lorentz force. For a 
"short" SNS junction (with a coherence length much 
longer than the nanowire) the AC Josephson dynam- 
ics of the device, induced by a weak DC driving volt- 
age, eV < (A£^(0))^/Ao, can be expressed in terms 
of a pair of Andreev states, which are periodically cre- 
ated and destroyed in the junction. Here, AE{(j)) is the 
distance between the Andreev levels, = (/)(t) is the 
voltage-driven phase difference over the junction and Aq 
is the order parameter in the superconductors. The cre- 
ation/destruction of the Andreev states takes place at 
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FIG. 1: (Color online) Schematic diagram of the system con- 
sidered. A suspended carbon nanotube of length L is coupled 
to two superconducting leads biased at a voltage V. Coupling 
between the Andreev states for the electronic degrees of free- 
dom and the mechanical vibrations of the nanowire is achieved 
through the external transverse magnetic field H which en- 
ables transition between the electronic branches through the 
emission/absorption of a vibrational quantum (see text). 



the edges of the superconducting gap, AE — 2Ao, when 
(j){t) = 2TTn (see Fig. [5]). Between creation and destruc- 
tion the Andreev levels evolve adiabatically within the 
energy gap of the superconductors, such that the inter- 
Andreev level spacing first shrinks and is later restored 
to its initial value 2Ao at a time when the Andreev states 
dissolve in the continuum quasiparticle spectrum (see 
Fig. [5]). The whole cycle then repeats itself again and 
again. 

Being thermally populated at the moment of their nu- 
cleation, the Andreev states will initially, as the inter- 
level spacing decreases, experience an effective cooling 
given that the thermal relaxation is not fast enough to 
follow the level displacement. Under such conditions, the 
Andreev states can accumulate energy from any exter- 
nal system (in particular from the nanomechanical vibra- 
tions) if interactions with such a system are established. 
The absorbed energy is then transferred into the elec- 
tronic quasiparticle continuum when the Andreev states 
are dissolved and the process of cooling of the exter- 
nal system is continued through the formation of new 
thermally populated Andreev states, which in their turn 
evolve with time and absorb energy. After each full pe- 
riod of evolution the extra energy absorbed in this process 
is removed from the SNS junction by the quasiparticle 
excitations. 

Besides being a new type of system where efficient 
cooling of nanomechanical-resonators can be achieved, 
this type of superconducting weak links makes it pos- 
sible to couple the nanomechanical system of the device 
to the macroscopically coherent superconducting states. 
It presents an exciting possibility to entangle the current 
carrying Andreev states associated with the SNS junction 
to the quantum states of nanomechanical vibration when 
the mechanical subsystem is sufficiently cooled. Manipu- 
lation of such nanoelectromechanical coherent states can 
be done both electronically, by means of driving- and gate 
voltages, or mechanically, by controlling the displacement 
of a movable part of the device. 
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2. MODEL HAMILTONIAN 

To discuss the cooling mechanism quantitatively we 
use the Hamiltonian, 



(1) 
(2) 
(3) 



In 111), Hei{t) is the Bogoliubov-de Gennes (BdG) Hamil- 
tonian for the electronic degrees of freedom. Here, 
A{x) = Ao0(|a;| — L/2) is the order parameter in the 
superconductive leads (here taken to be identical), fi is 
the chemical potential and the phase difference over 
the junction, depends on the bias voltage V according to 
the Josephson relation = 2eV/h. The second term in 
the Hamiltonian describes the oscillating nanowire which 
is modeled as a quantum mechanical harmonic oscillator 
with only the fundamental mode excited; [b] is an op- 
erator that creates [annihilates] a quantum of vibration 
and u is the frequency of the fundamental mode. In ([T]), 
Hint is the interaction Hamiltonian which gives the cou- 
pling between the electronic and the mechanical degrees 
of freedom through the magnetic field. 

At the superconductor-nanowire interface two elec- 
tronic scattering processes are possible, viz. normal re- 
flection and Andreev reflection. These give rise to two 
localized states — Andreev states — described by the 
wavefunctions which satisfy Heiip±ix,t) — E±'ip±{x,t) 
where 'i{}±{x,t) are two-state spinors in Nambu space, 
whose eigenenergies are given by the expression E± (t) — 
±Ao(l-i:'sin2((/)(t)/2))i/2^here D = 1 - i? is the nor- 



mal junction transparency TsUTqI . These Andreev states 



are responsible for carrying the Josephson supercurrent 
through the (normal) nanowire. Under the condition that 
these states evolve slowly in time {h(j) ^ Aq), they de- 
scribe a two-level system for the electronic degrees of free- 
dom. As the Andreev states carry current through the 
oscillating nanowire, they will couple the electronic de- 
grees of freedom to the mechanical degrees of freedom 
through the Lorentz force. It is thus more natural to 
describe the system in terms of the temporally evolving 
Andreev states, to which purpose we rewrite the Hamil- 
tonian as 



H{t) = <^{<P{t))a, + AoVi?sin(0(t)/2)a,+ 



hajb^b- 



2e dS{(t){t)) 



Here, S'{(p{t)) is the energy of the Andreev states for 
the completely transparent (i? — 0) junction and ai 
{i = x,y, z) are the Pauli matrices. The last term in (|4]) 
describes the electromechanical coupling of the current- 
carrying Andreev states to the motion of the nanowire 
through the Lorentz force. Here, e is the electronic 
charge, L is the length of the wire, y = yo{b'' + b) — 



Ao 





'■■<. 




o 








v" " ' ■ 


!/2{?i)P|-.ra) 


o ^ 


/ \d 


[n)\^\-.n)/ 


to Ty 



FIG. 2: Time evolution of the Andreev states (full lines) over 
one period Tv = TxhjeV . The state of the total system, |±, n), 
depends on the population of the two electronic branches, |±), 
corresponding to the upper and lower branch respectively, as 
well as on the quantum state of the oscillating nanowire, \ri). 
Due to the large separation in energies between Ao and fttj, 
transitions between the electronic branches is only possible in 
the small resonance-window to — (5t < t < to + i5t (modeled 
through the scattering matrix S) where the electronic state 
of the system can change through the emission/absorption of 
one vibrational quantum. In the above, the probability for the 
state initially in | — , n) to scatter into the state |-|-, n — 1) after 
passing through the resonance depends on the state of the 
oscillator through the coefficient |i/2(»i)P (see text). After one 
period the partially filled Andreev levels join the continuum, 
a process which is here represented by dashed arrows, and the 
electronic states are reset (filled and empty circles). 



where yo is the zero-point oscillation amplitude — is the 
deflection operator in the transverse y-direction, H is the 
magnetic field and 2e[h{dS / d(t))<Tz is the current opera- 
tor in the nanowire l20| . 

From Fig. [5] one can see that the cooling mechanism 
described is maximally efficient if the strongest coupling 
of the mechanical subsystem to the Andreev states occur 
during the time interval when the Andreev levels are clos- 
est to the chemical potential, i.e., when the population 
of the Andreev levels differ most from what it would be 
in thermodynamic equilibrium. The coupling is resonant 
during this time if hiu = Ttvai{/S.E{4>)) = 2Ao\/i? (where 
R = 1 — D \s the normal-state reflection probability), 
which requires the weak link to have a high normal-state 
transparency D for electrons [2l|. This is because, for 
the system considered, the energy-scales of the super- 
conductive order parameter and the mechanical vibra- 
tions, tujj, are very different. Also, one notes that when 
(4) = TT, the Andreev states ip±{4> — tt), with energies 
E±{4> = tt) = ±Aoi?^/^, are symmetric/antisymmetric 
superpositions of states carrying current in opposite di- 
rections. Transitions between the Andreev states — in- 
duced by the nanowire as it vibrates in a transverse mag- 
netic field — are therefore by far most probable when 
0(i = ^o) TT. 

In the adiabatic regime, i.e. when eV < eVc = 4i?Ao, 
it is convenient to switch to the basis set formed by the 
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Andreev states il^±{4>{t)^^. In this basis the Hamiltonian 
reads, 

Heffit) = {iJ±\H{tM±) = E{t% + 

fta;6t6 + Ao$(6^ + fe)f^. (5) 

Here, the Pauh matrices fi span the space formed by the 
states -(Alt (</>)■ In the last term of (O, $ = 2LHTryo/^o 
is the dimensionless magnetic flux threading the area 
swept by the nanowire as it vibrates in the ground state; 
$0 = h/2e is the magnetic flux quantum. Below we 
will consider the resonant situation when hu) = 2Aoi?^^^, 
which is the optimum condition for the proposed cooling 
mechanism as outlined above^^. 

The Hamiltonian ([5]) describes the evolution of the An- 
dreev states. Coupling between these is achieved through 
the off-diagonal terms in ^ . These are proportional to $ 
and become relevant whenever the energy separation be- 
tween the electronic branches matches the energy scale 
set by the oscillating nanowire. At such times, transi- 
tions between the Andreev states can occur through the 
absorption/emission of one vibrational quantum. This is 



shown in Fig. [5] for the case when only the lowest elec- 
tronic branch is initially populated. 



3. MAGNETIC FIELD INDUCED ANDREEV 
STATE SCATTERING 

Assuming the conditions for adiabaticity outlined 
above fulfilled, we now evaluate the probability of me- 
chanically induced Andreev level transitions as a function 
of the coupling strength. Due to the large separation of 
energy scales, 2Ao » huj, it is not possible to evalu- 
ate such transitions using the Landau-Zener formalism 
20||2^ as this assumes a linear expansion of E{t) around 
the resonance point. Rather, we consider transitions only 
to occur in the vicinity of to — Tv/2 ~ TTh/2eV and eval- 
uate ^ using a parabolic expansion of E{t). From ([SJ, 
the dimensionless time-dependent Schrodinger equation 
for the probability amplitudes c±^„(t) of finding the state 
of the system in the upper/lower electronic branch with 
the oscillator in the state n reads. 



i9rC-,„(r) 



(6a) 
(6b) 



(t - to) 



1/3 



£Ao 
hoj [v 



2/3 



D{hjjjf 



Xt 



2/3 



T:huj 
2eV 



Equation ([6]) describes the adiabatic evolution of the 
electronic Andreev states. As theses states pass through 
the resonance they couple to the mechanical subsystem, 
enabling transitions between them through the emis- 
sion/absorption of a quantum of vibration as indicated 
by the terms proportional to the dimensionless electron- 
vibron coupling constant F in ([6]). Considering for now 
the situation when the system initially starts in the lower 
electronic branch (I6ap we analyze the probability that the 
system, after passing through t^, scatters into the upper 
electronic branch as a function of the coupling strength. 
As can be seen from ([5a]) . two scattering processes are 
possible; absorption of a mechanical vibron (proportional 
to yfn) or emission of a vibron (proportional to \/n — 1). 
Note however, that these two processes differ by an expo- 
nential phase factor associated with the emission process. 

Considering that the coupling between the electronic 
and mechanical degrees of freedom is very weak (see be- 
low) we evaluate the probability for the electronic sub- 
system to scatter into the upper electronic branch per- 
turbatively in the small parameter F. From this analy- 
sis we find that scattering is much more likely to occur 



through the absorption process than through the emis- 
sion process, whose importance is reduced by the expo- 
nential phase factor mentioned above. A numerical anal- 
ysis shows that the probability of magnetic-field induced 
Andreev state scattering through the emission channel is 
only ^ 3% of the probability of scattering through the 
absorption channel. For the case when the system is ini- 
tially in the upper electronic state ((6b| the same analysis 
shows that the emission process (proportional to ^/n) is 
much more likely than the absorption process (propor- 
tional to ^/n + 1). Thus, we can safely treat the absorp- 
tion process (when the electronic system is initially in 
the lower electronic branch) as the most dominant pro- 
cess, i.e. we may apply the rotating wave approxima- 
tion (RWA) to ([B]). Note, however, that without thermal 
damping of the nanowire vibrations caused by the exter- 
nal environment (see below) , the rate of emission would 
ultimately limit the efficiency of our cooling process as 
it would prohibit complete ground state cooling. The 
details of such an analysis is presented in Appendix |Bl 
where we show that by not using the rotating wave ap- 
proximation the lowest theoretically achievable level of 
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cooling of the nanowire corresponds to a population of 
the vibrational mode of (n) ~ 0.03. 

Within the rotating wave approximation we proceed 
to evaluate the scattering probability as a function of the 
coupling strength in more detail, only keeping the terms 
proportional to ^Jn in (jS]), 



idrCnir) (r^T^ + T^/nr^) c„(t) 
c„(t) = (c+^„_i(t),c_^„(t))'^ . 



(7) 



The electronic scattering mechanism considered here is 
similar to the Landau-Zener formalism in that the prob- 
ability of scattering depends on the ratio between the 
coupling strength and the rate of evolution through the 
interaction region, off-diagonal terms in However, 
unlike in the Landau-Zener approach, where in the cor- 
responding equation the dependence on r is only linear, 
equation ^ cannot be solved analytically. Instead we 
have performed numerical simulations in order to find 
the probability of scattering as a function of the cou- 
pling strength. From this analysis we find that unlike the 
Landau-Zener formula, which gives a unitary transition 
probability in the strong coupling limit, the parabolic 
dependence on r in (O dictates that the probability of 
transition in the limit of strong coupling (or infinitely 
slow transition through the interaction region) is only 
0.5. Thus, for very low applied bias voltages we should 
expect a probability of transition between the Andreev 
states of 0.5 each time the system passes through the 
resonance. 

In the adiabatic regime, V <Vc the coupling for realis- 
tic experimental parameters is, however, relatively weak; 
F = Ao$/(fia;) ^ 1, hence we can analyze Q treating F 
as a perturbation to the solution |c_,„(— 5t)P — 1, where 
the system starts in the lower electronic branch at time 
-St = -St{^/hy/^. Here, St ~ {huj<^>/Aoy^^h/ieV) < 
ttH/ (eV) is the characteristic time scale of the electrome- 
chanical interactions. Solving ^ under this condition 
for the probability |c+,„_-i((5t)P of the system to be in 
the upper branch after passing through the resonance we 
find |c+^„_i((5r)p ~ mrV^ which is consistent with our 
numerical solutions. Thus, for weak coupling between the 
electronic and mechanical degrees of freedom, the proba- 
bility for a transition between the Andreev states to have 
occurred after passing through the resonance (with the 
corresponding absorption of a mechanical vibron) scales 
linearly with the initial population of the mechanical res- 
onator. Furthermore, as can be seen from the symmetry 
of ([7]) , this expression also gives the probability of finding 
the system in the state \ — ,n) given that it was nucleated 
in |+, n - 1). 



4. GROUND STATE COOLING OF 
OSCILLATING NANOWIRE 

In order to model the full evolution of the coupled elec- 
tromechanical system we evaluate the total density ma- 



trix p of the system over one period. Under the assump- 
tion that the thermal energy at the ambient temperature, 
T, is small compared to the initial separation between 
the electronic branches, 2Ao/fcsT > 1, the system will, 
at the start of the period, be found in the lower elec- 
tronic branch with a distribution, , of the mechanical 
excitations, 



i,7=±ri=0 n=0 ^ ^ 



\n){n\ . 



During the adiabatic evolution (no coupling between 
branches) the system interacts with the external heat 
bath and the rate of change of the density matrix is given 

by, 



dm 

dt 



H,fAt),m\+jmt)): (8) 



where the collision integral 

ne (bPp + pSSt - 2pph) , (9) 



models the interaction of the mechanical subsystem with 
the environment. Here, Ub = {exp{hw / ksT) — and 
7 = w/Q is the thermal damping rate of the vibrational 
modes with Q the quality factor. Considering the cou- 
pling to the environment to be small {Q ^10^ [1^) we 
can solve for the evolution of the density matrix by treat- 
ing the interaction of the system with the external envi- 
ronment as a perturbation. 

The adiabatic evolution of the electronic system en- 
sures that the coupling between the Andreev states is 
virtually zero. Only in the small region around tg do 
the two branches interact (through the coupling to the 
mechanical degrees of freedom) , which is conveniently ac- 
counted for by describing the evolution of p through the 
resonance by the unitary scattering matrix S, 



p{to + St) = Sp{to - 6t)S^ 



Ki{h) 



\K^{n)\ 



I'ljn) ^ 
2 I I,, f„\\2 



(10) 

(11) 



1 



Here, fi — 6^6 is the vibron number operator and the sub- 
scripts 1, 2 refer to the top/bottom Andreev state respec- 
tively. The coefficients Ki(n) [i'i(fi)] are respectively the 
probability amplitude for the system to stay [scatter] out 
of the initial Andreev state i, both of which depend on the 
state of the oscillator as outlined above. As such, |i^2(?^)P 
is the probability of the system, initially in the state n 
in the lower electronic branch, to scatter into the upper 
electronic branch through the absorption of a vibron. It 
thus corresponds to |c+,„_i((5r)p from ^ with the initial 
state | — ,".). With this we have |i^2(f^)P — ttoF^ and also 
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of 0.15 



FIG. 3: (Color online) Evolution of the distribution of the 
mechanical modes, P„, as a function of the quantum state 
n for different number of periods {Tv ~ nh/eV ~20ns). 



Initially P„ is thermally distributed Pn (x e 



-hujn/ k 



with 



=10" 



eV, Ao = Whw, 

2 „„ 



kgT = 5huj. Here, V =10"' V, hw 
yo =20pm, L =100nm and H =1 T. Inset shows \T^2{n) 
a function of n for the same parameters. 



|i^2('T-)P = Wiin ~ 1)P which follows from the symmetry 
of (0). 

After one period the partially filled Andreev branches 
merge with the superconducting continuum. Here, non- 
adiabatic interactions release the energy of the charge- 
carrying quasiparticles into the continuum and the initial 
conditions for the Andreev level population is reset p^ . 
i.e. the electronic system returns to the lower branch 
(see Fig. [2])- The mechanical system, which forms the 
major interest of this paper, will however not return to 
the initial distribution after one period. Thus, we may 
find the density matrix for the mechanical system after 
one period, pf„gc/i' ^"7 tracing out the electronic degrees 
of freedom of the total density matrix. 



n=0 



-J - 



(12) 



Equation (|T2|) describes the evolution of the mechani- 
cal density matrix over one period. Without coupling to 
the environment, 7 = 0, this corresponds to the decay 
of the mechanical subsystem as for each period there is 
a probability |j/2(?t-)P that the state \~,n) scatters into 
the state \ + ,n ~ 1). As the opposite process is forbid- 
den if initially only the lower electronic branch is pop- 
ulated, the mechanical subsystem would thus approach 
the vibrational ground state at a rate which depends on 
the strength of the coupling. As noted above, complete 
ground state cooling would however not be possible due 
to the small but finite probability of inter- Andreev level 



scattering through the emission channel. 

When including effects of the environment, the me- 
chanical subsystem does not necessarily decay to the 
ground state since the external damping drives the sys- 
tem towards thermal equilibrium. In Fig. [3] we plot the 
evolution of the distribution of the vibrational states, 
Pn, as a function of the quantum numbers n over many 
periods for realistic experimental parameters. As can 
clearly be seen, the above discussed scattering process 
acts to lower the initially thermally distributed popu- 
lation of the mechanical system, effectively cooling the 
nanowire down to a final population of the vibronic states 
of (n) = Y^n^^n ^0.1. Note, however, that complete 
ground state cooling is not achieved due to the tempera- 
ture of the external environment. 

To evaluate the efficiency of the discussed cooling 
mechanism we analyze the stationary solution of (I12p . i.e. 
when the mechanical subsystem does not change over one 
period, perturbatively in the small parameter 7. From 
this analysis we find that to zeroth order in 7 the sta- 
tionary solution gives \i'2(n)\'^Pn = which corresponds 
to the distribution of the mechanical modes decaying to 
the ground state (|i^2(0)P = 0). The corrections to this 
solution due to the external heat bath is found from the 
higher order terms in the pertubative expansion where 
the lowest order correction to the state n scales as 7". 
We thus find that to linear order in 7 the correction to 
the distribution of the mechanical modes in the station- 
ary regime, P^'"*, is given by 



ps 



SPn = (-1) 



TTftns 

2eV " 

Pn — PqSu.O 

„+i (2-|^2(l)P) 

l^2(l)P 



0(7') 



71 = 0,1, 



(13) 



With this, we find that the population of the first excited 
state in the stationary regime can be expressed as. 



psi 



1/3 



0.1, 



(14) 



which is in accordance with the results shown in Fig. [31 
In this analysis, the strength of the electromechani- 
cal coupling, r, plays two roles. First, it dictates the 
rate of the cooling process as |i/2(n)p cx F^. Secondly, 
stronger coupling ensures a lower final distribution of the 
mechanical modes; in effect the ratio 7/r^ dictates how 
many terms in the perturbation expansion (1131) need to 
be considered in order to accurately describe the station- 
ary distribution of the mechanical modes. Thus, for the 
parameters considered, we find that we need to include 
higher order corrections in 7 in the expansion (jl3p to fully 
describe the stationary distribution. This is confirmed 
in Fig. [3] which shows a finite probability of not only 
the ground state and first excited state being populated 
in the stationary regime. Nevertheless, the expectation 
value of the population of the mechanical subsystem in 
Fig. [3] is (n) ~ 0.1, which is well in the quantum regime. 
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5. DC CURRENT MEASUREMENT AS A 
PROBE FOR THE STATIONARY VIBRONIC 
POPULATION 

To probe the stationary distribution of the mechanical 
subsystem we suggest studying the DC current through 
the weak hnk as this directly measures the vibron popu- 
lation as will be shown below. By changing the quality 
factor of the system, for example by varying the external 
pressure, we find that the stationary vibron population 
can be directly measured since it scales linearly with the 
DC current. To show this we evaluate the DC current, 
on resonance, induced by the inter- Andreev level scat- 
tering over one period when the mechanical subsystem 
has been driven into the stationary regime. This current 
comes about as the two Andreev levels carry current in 
opposite directions. As such, any scattering event which 
populates the upper Andreev level will result in a net 
charge transfer through the system. Evaluating the total 
current transfer over a full period we have. 



2e 2 



.^0 



JTv/2 



JTv/2 



where, P is the total probability for the system to be 
scattered from the lower to the upper Andreev branch 
during one period of evolution. Since n vibration quanta 
are excited with probability P,f*°*, and since in this case 
the scattering probability is |i^2("-)P, the total scattering 
probability P is readily obtained, and the expression for 
Idc simplifies to (note that the time integration is trivial 
since d/dcj) = {dt/d(t))d/dt = {Tj2n)d/dt) 



4eAo A hu 



71=0 



h \ R 



(16) 



which is plotted in Fig. 2] for the same experimental pa- 
rameters as in Fig. |3l 

Equation (fT6|) shows that the DC current over the junc- 
tion scales linearly with the average vibron population 
in the stationary regime. Here, we note that the cor- 
rect expression for the current should also take into ac- 
count corrections from the environment as discussed in 
Appendix [C] Numerical simulations show, however, that 
the inclusion of these terms does not change the calcu- 
lated DC current as compared to that shown in Fig. |4l 

From Fig. |4] we note that in the limit of small quality 
factor, i.e., high damping, the DC current saturates to 
a constant values. This can be understood by the fact 
that in this limit the stationary distribution of the vi- 
brational modes is simply given by the thermal distribu- 
tion P^*"* c>c exp{—nhuj/kBT), i.e. the interactions with 
the environment are so strong that they always drive the 




10" 



10' 



FIG. 4: DC current as a function of the inverse quality fac- 
tor when the mechanical subsystem has been driven into the 
stationary regime. As can be seen, in the limit of very high 
quality factor, corresponding to complete ground state cool- 
ing of the oscillating nanowire, the DC current goes zero as 
no inter-Andreev level scattering is possible in this regime 
(|z^2(0)P = 0). In the opposite limit, Q < 10, the DC current 
approaches a constant value which depends on the external 
temperature, I^^^'" cc ub- Here, the system parameters are 
(15) the same as in Fig. [3l 



mechanical system into thermal equilibrium, no matter 
how the inter-Andreev level scattering changes this dis- 
tribution. As such, I DC oc Ub in this regime. In the 
opposite regime, Q — >■ oo, the DC current goes to zero 
(within the rotating wave approximation) as expected as 
this regime corresponds to complete ground state cooling 
of the mechanical subsystem. If this can be achieved, the 
probability of Andreev level scattering also goes to zero, 
|i^2(0)P = 0, hence the electronic subsystem stays in the 
lower Andreev branch throughout and the loc = 0. 

From the dependency of Idc on l/Q in Fig. [3] we sug- 
gest a simple probe by which one could measure the av- 
erage vibron population once the mechanical subsystem 
has been driven into the stationary regime. We propose 
to change the quality factor, for example by increasing 
the external pressure of the gas in which the nanowire vi- 
brates, and measure the DC current in the limit of large 
damping. As Idc scales with the thermal equilibrium 
phonon population, ng, in this regime a further mea- 
surement in the low damping regime should thus give a 
direct measurement of the average vibron population. 



(n) 



Idc 



jtherm 
^DC 



(17) 



where I^^^"^ oc ub is the DC current in the limit of high 



damping. 
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6. CONCLUSIONS 

To conclude, we have shown that quantum mechanical 
cooling of a nanomechanical resonator acting as a weak 
link between two superconductors is possible. In par- 
ticular, we have considered the example of a suspended 
carbon nanotube where the unique combination of high 
resonance frequencies and mechanical quality factors to- 
gether with high electronic transparency lead to suffi- 
ciently strong coupling for efhcient cooling to be possible. 
Using realistic experimental parameters we have shown 
that for a short suspended nanotube, stationary distribu- 
tion of the population of the vibrational modes can be as 
low as (n) =0.1 which is truly in the quantum mechanical 
regime. Furthermore, the suggested mechanism does not 
rely on any external electromagnetic fields to stimulate 
the cooling process. Rather, the proposed system should 
act as a self-cooling device, given the right experimental 
parameters, as the over-cooled Andreev states can readily 
absorb energy from the mechanical subsystem given that 
sufficient coupling between the two can be achieved. The 
corresponding energy uptake of the electronic subsystem 
is later released into the quasiparticle continuum, leading 
to an effective cooling of the nanomechanical resonator. 

Furthermore, it has been shown that by measuring the 
DC current as a function of the quality factor, the sta- 
tionary distribution of the population of the vibrational 
modes can be directly probed. Thus, the suggested sys- 
tem not only has the capacity of self-cooling, but it also 
presents a direct, internal probe by which to measure the 
efficiency of the cooling mechanism. 
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Appendix A: Electrical field induced cooling of 
suspended nanowire 

As mentioned above the required coupling for effec- 
tive cooling of our suspended vibrating nanowire could 
also be achieved using a transverse electric field. Un- 
der such situations, the non-superconducting suspended 
wire supplies a coupling of the Josephson current flow 
to the nanomechanical bending vibrations through the 
displacement of the wire with the electromechanical cou- 
pling depending on the electronic backscattering induced 
by the electrical field. With this setup, the interaction 
term in the Hamiltonian ([1} reads, 

H^nt = eSy^iP + b)a, , (Al) 



where S' is the transverse electric field. Proceeding as 
above, we project out the Andreev states -^i, which gives 
the effective Hamiltonian, 

H,f}{t) = {ip±\Hm±) = AS(<)f, + nwb'^b+ 

S(St+6)(f+U+_+h.c.) (A2) 
l^E{t) ^ E{t) + \{y^ +h)V++ E^eS'yo. 

In the above, Vij = {Tpi\az\ipj) , t_|_ — {t^ + ify)/2 and 

I oc \/R [12]. To evaluate this further we perform 

the unitary transformation UHef f{t)U\ 

U = exp(a(6t - b)fz) a = EV++/huj < 1 , (A3) 

after which the Hamiltonian reads, 

Jifeffit) ^ E{t)Tz + hujb^b+ 

S(6t+S)(U+_f++h.c.) + 0(a2). (A4) 

From here on, the analysis follows the previous 
one with the coupling constant now being, F — 
\/hLij){Vc/V)^^'^ ■ For realistic experimental pa- 
rameters, our analysis shows that the coupling constant is 
again much smaller than unity. Hence we can develop the 
same perturbative analysis for the evolution of the den- 
sity matrix. Using the same physical parameters as above 
for the dimensions of the nanowire, the external tempera- 
ture and the superconductive gap with =4x10^ Vm^^ 
we find a stationary occupancy of the vibrational modes 
of the oscillating nanowire of (n) ~ 0.25. Thus, this anal- 
ysis shows that an alternative approach to that followed 
in the main part of the paper would be to use an elec- 
tric field to provide the electromechanical coupling. Both 
modes of operation give an efficient cooling of the vibrat- 
ing nanowire. However, our analysis shows that the best 
mode of operation for ground state cooling should be to 
consider using a magnetic field. This mode also has the 
added benefit that it makes it possible to tune the su- 
perconductive order parameter with the magnetic field 
to achieve resonant conditions. 



Appendix B: Emission limited cooling of the 
nanowire vibrations 

The analysis presented in the main part of the pa- 
per, performed within the rotating wave approximation 
(RWA), showed that the suggested cooling mechanism 
has the capacity of effectively reducing the excitations of 
the nanowire vibrational modes well into the quantum 
regime. Within the RWA, complete ground state cooling 
is, in theory, achievable as a large enough quality fac- 
tor limits the influence of the environment (heating) on 
the mechanical sub-system. Thus, in the limit Q — >■ oo, 
thermal heating is negligible and the inter- Andreev level 
scattering mechanism reduces the vibronic population to 
the ground state with the corresponding decay of Idc- 
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In the rotating wave approximation we ignore the pos- 
sibihty of electronic scattering from the lower to the up- 
per branch through the emission channel. This process 
can occur if the time spent on resonance, 6t, is short in 
which case the uncertainty principle dictates that scat- 
tering can occur not only through the absorption but also 
through the non-energy conserving emission channel. By 
considering also this channel for inter- Andreev level tran- 
sition, complete ground state cooling is not achievable 
(even in the limit Q — >■ oo) as the emission process will 
continuously heat the oscillating nanowire. To analyze 
this we consider ^ in more detail. 

Without the benefit of the RWA approximation we 
need to consider all four off-diagonal terms in ^ in order 
to fully describe the dynamics of the scattering mecha- 
nism. From the form of (I6al) one can see that e.g. the 
term proportional to \/n — 1 corresponds to the scatter- 
ing from the lower electronic branch in the state n — 2 
to the upper electronic branch in the state n — 1 through 
the emission of a vibron. Thus, this rate of scattering 
depends on the probability of the occupancy of the har- 
monic oscillator state n — 2 and is proportional to n — 1 
as compared to the absorption channel which, as shown 
above, depends on the probability of the population of 
the state n and scales linearly with n. To take these ef- 
fects into account we rewrite the unitary scattering ma- 
trix nil), 



K[{h) 



" '-1 \/rL 



/n+T 




(Bl) 



In (IB1[) I K2('T-)P, 1 1^2(^1)1^ have the same meaning as 
before, i.e. they are respectively the probability to stay or 
scatter (through the absorption channel) from the lower 
electronic branch with oscillator in the state n. Similarly 
|X2(f^)P gives the probability that the system scatters 
from the lower to the upper electronic branch through 
the emission channel (the subscripts i = 1 correspond to 
scattering from the upper to the lower electronic branch). 
By construction, these rates are, 

|^2(n)P = |^i(n-l)p 
\Xi{n)?^\X2{n-l)? 



|X2(n-2)|2 „ 



\v2{n)[^ 



n — 1 



0.03, 



where |i^2('T')P = mrT'^ and the coefficient of n/{n — 1) 
in the last line was found from numerical analysis as dis- 
cussed above. 

To analyze the limiting value of the stationary distribu- 
tion of the vibrational modes in the non-RWA formalism 
we consider the evolution of the mechanical sub-system 
with Q 00. Proceeding as above we evaluate with 
3 = 82 and 7 = 0. From this we find that the distribu- 



tion of the mechanical modes after one period reads. 

Pi ^ PniiW2{n+i)\'+p:r\^'2{n)\'+ 



C-ilX2(n-l)|^ 



(B2) 



Equation (|B2p dictates that complete ground state cool- 
ing of the mechanical sub-system is not achievable as the 
term proportional to |x2(?^ — 1)P acts to shift the popu- 
lation to higher vibrational modes and will thus compete 
with the term proportional to 1^2 (n + 1) which acts to 
lower these. Considering however, that the rate of ab- 
sorption is much greater than the rate of emission we 
find that the stationary solution of (|B2p corresponds to 
a population of the vibrational modes corresponding to 
(n) — 0.03 which is in accordance with the perturbative 
solution, 



\X2[ 



0.03. 



(B3) 



Thus, we find that by including not only the absorption 
but also the emission channel for inter- Andreev scattering 
the best theoretically achievable level of cooling is given 
by the ratio of these two rates. The corresponding cur- 
rent with this population of the vibrational modes would 
be Idc —l.l pA. 



Appendix C: Environmental correction to the DC 
current 



As discussed in the main text the DC current through 
the nanowire scales with the average vibron popula- 
tion, {n). In this is shown for the case when 
P = J2nW2{n)\^P^*°-'^ . We note here that this expres- 
sion is not completely valid due to the influence of the 
external environment. To understand why this is the 
case we have to consider that the stationary distribution 
pstat gjygg vibronic population which does not change 
over one full period. Over this interval of time, the me- 
chanical sub-system changes in two ways; by temporal 
evolution and through magnetic-field induced scattering. 
The value of P^*"^* depends on both of these processes. 
This means that in order to accurately evaluate the total 
probability of scattering of the mechanical sub-system we 
first have to consider the effects of the environment on 
the evolution of P^*°* between < t < Tv/2 and calcu- 
late the total probability of scattering from this quantity. 
Performing this analysis we find that the accurate prob- 
ability of scattering reads. 



P = vrF^ 



(CI) 



which is again linearly dependent on (n) but also includes 
a correction term from the environment. 

A numerical analysis of the DC current as a function 
of Q has been performed using both (IT5| and (jCll) . with 
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little or no discrepancy between the two results. To un- 
derstand why the corrections in (|Uip does not change 
I DC we analyze separately the limit of low/high damp- 
ing which both give P = 7rr^(n). In the limit of high 
quality factor, the corrections to i-"^*"^* from the thermal 
interactions are very small cx l/Q, hence they do not in- 
fluence the stationary mechanical distribution over the 
short time-span Ty/2 to any large extent. On the other 



hand, in the limit of low quality factor the stationary dis- 
tribution of the mechanical modes is given by the thermal 
distribution, hence the correction to the total scattering 
probability in (ICip are identically zero. For intermediate 
values of Q these competing processes cancel out and the 
current over the junction depends only on the average 
vibron population in the stationary regime as given by 
equation pB)). 
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tude of the order parameter Ao through the magnetic field 
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